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The compositional dependence of thermal expansion behaviour in 19 different perovskite-like metal-organic frameworks (MOFs) 
of composition [A^][M^^(HCOO)3] (A = alkylammonium cation; M = octahedrally-coordinated divalent metal) is studied using 
variable-temperature X-ray powder diffraction measurements. While all systems show essentially the same type of thermome¬ 
chanical response—irrespective of their particular structural details—the magnitude of this response is shown to be a function of 
A^ and cation radii, as well as the molecular anisotropy of A^. Flexibility is maximised for large and small A^ while the 
shape of A^ has implications for the direction of framework hingeing. 


1 Introduction 

The discovery that many metal-organic frameworks (MOFs) 
respond mechanically to external stimuli in extreme and coun¬ 
terintuitive ways {e.g. ne^tive thermal expansion (NTE),^ 
negative conipres sibilityj^ breathing transitions,® and 
amorphisation®) has focussed attention within the field on 
establishing composition-property relationships. The hope— 
as in much of MOF science—is to develop design strategies 
that will allow the targeted synthe sis of MOFs with specific 
and optimised physical properties.Given the enormous 
structural and compositional diversity of this family, the task 
of establishing design rules has focussed on two separate ar¬ 
eas: namely, establishing in turn the roles of network topol¬ 
ogy/geometry and network composition in governing elastic 
response. 

With regard to the former, it is now well established that 
certain network connectivities intrinsically favour counterin¬ 
tuitive mechanical responses Perhaps the best known 

example is the wine-rack topology, which is associated 
with anomalous elastic behaviour across a large variety of 
chemically-distinct MOF familiesFor a given network 
topology, changes in geometry {e.g., network angles) can in- 
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vert the type of a nisotr opy but the basic mechanical response 
remains the same.^EH 

With regard to network composition, variation of the 
chemical makeup of MOFs with the same topology will 
usually infiu ence only the magnitude of mechanical re¬ 
sponse. ®E 3 E 11 For instance, nanoindentation studies of the 
perovskite-structured [(CH3)2NH2][M(HCOO)3] family (M 
= Mn, Co, Ni, Zn) showed a strong relationship between 
ligand field stabilisation energy of the transition metal di¬ 
cation and the mechanical stiffness of the corresponding 
MOF.ISlln the related systems [C(NH2)3][Mn(HCOO)3] and 
[(CH2)3NH2][Mn(HCOO)3], the hydrogen-bonding strength 
of the extra-framework alkylammonium cation was found to 
direct the stiffness and fiexibility of the two structures.in 
all cases, the fundamental mechanism of elastic response is 
unchanged by chemical substitution, with the degree of fiexi¬ 
bility scaling inversely with strength of interaction .El Precisely 
the same conclusion has been reached in similar studies of in¬ 
organic frameworks. For example, an investigation of the ther¬ 
mal expansion response of Prussian Blue analogues and re¬ 
lated lanthanide hexacyanocobaltates showed that cation size 
correlates with the magnitude of NTE.E 3 E 2 I Likewise, cation 
size also plays an important role in the magnitudes of thermal 
expansion behaviour in the family of frameworks related to 
NaZr2(P04)3.E3 Perhaps the only exceptions to these general 
rules are in instances where chemical substitution alters the 
nature o f the d ominant chemical interaction within a particu- 
lar MOF.iniSl 

One particularly important way in which dense MOFs can 
differ from “conventional” ceramic frameworks is in their ca¬ 
pacity to accommodate extra-framework ions that are molec¬ 
ular rather than monatomic in nature. The added complexity 
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[C(NH2)3] [Co(HCOO) 3] BaTi 03 


Fig. 1 ABX3 frameworks of [C(NH2)3][Co(HCOO)3] and BaTi03. 
The hydrogen bonding in [C(NH2)3][Co(HCOO)3] is shown with 
dotted red lines. The Co^^ and Ti^+ coordination environments are 
represented by polyhedra. 


of cation asphericity is conceptually related to the symmetry- 
lowering effect of second-order Jahn-Teller instabilities and 
may play a key role in the lattice dynamics of high-profile 
hybrid frameworks such as [CH3NH3]Pbl3.ISl Yet, to the best 
of our knowledge, there are no systematic studies of the rela¬ 
tionship between counterion shape and mechanical response 
in MOF-type systems. 

In this paper, we use thermal expansion measurements 
to explore the relationship between framework fiexibility 
and chemical composition across 19 MOFs drawn from 
the widely-studied family of alkylammonium transition-metal 
(“ABX3”) formates. These compounds have the same general 
formula [A^][M^^(HCOO)3] and are structurally related to the 
perovskites [Fig.[T|. The larger M.. .M separation in formates 
relative to perovskites (^6 A 4 A) allows incorporation of 
molecular cations within the the framework cavities. As in the 
perovskites there is scope for substitution on both the “twelve- 
coordinate” A-site and on the transition-metal B-site. While 
the cation charges on each of these sites are fixed at 1 -F (A) and 
2 + (B) in the formates (unlike perovskites), there is now an ad¬ 
ditional degree of freedom in terms of the molecular shape at 
the A-site. 

Their relative structural simplicity makes these compounds 
ideal candidates for a composition/property study such as 
ours. But the broader family is also of significance from a 
functional materials viewpoint. ABX3 formates exhibit a va¬ 
riety of useful properties, including ferroelectricitY®Ell fer- 
roelasticity,!^ and even multiferroic behaviour.®® Compo¬ 
sition/property studies of the ferroelectric response in the A^ = 
NH4 or (CH3)2NH2^ members points to the different possi¬ 
ble effects of substit ution on the alkylammonium site and the 
transition-metal site.^^EE] Substitution at the A-site allows 
tuning of the ferroelectric polarisation,® whereas variation in 
tra nsitio n-metal affects the ferroelectric transition temperature 
j | 38 | 47 | context, framework fiexibility may have im¬ 


portant implications for A-site orientational ordering and/or 
mechanisms of accommodating the large strains induced dur¬ 
ing ferroelectric/paraelectric switching. So this study has 
the additional relevance beyond establishing composition- 
property relationships of exploring the role of fiexibility in 
ferroelectric MOFs. 

Our paper is arranged as follows. We begin by summaris¬ 
ing the synthesis and characterisation techniques used in our 
study. In the results section that follows, we report the ther¬ 
mal expansion properties of the various ABX3 formates we 
study. Using the mechanical building unit (XBU) abstraction 
developed in Ref. [H we reduce the experimental lattice pa¬ 
rameter data we measure to two characteristic values for each 
system: namely, the expansivities of the framework struts and 
of the intra-framework angles. By comparing the magnitudes 
of these two parameters as a function of framework composi¬ 
tion, we are able to establish composition/fiexibility relation¬ 
ships for this family. 


2 Experimental Methods 

2.1 Sample preparation 

The 19 structures we investigated share the composition 
[A’][M“(HC00)3], where = CH3NH3, (CH3)2NH2, 
CH3CH2NH3, (CH2)3NH2, C(NH2)3 and M" = Mg, 

Mn, Fe, Co, Ni, Cu, Zn, Cd [Fig. |^. For simplic- 
ity, the [CH3NH3][M(HCOO)3] structures will be re¬ 
ferred to as MeNH3-M, where M is the metal cation. 
Likewise we refer to [CH3CH2NH3][M(HCOO)3] as 
EtNH3-M, [(CH3)2NH2][M(HC00)3] as Me2NH2-M, 
[(CH2)3NH2][M(HC00)3] as Aze-M (i.e., Aze = azeti- 
dinium), and [C(NH2)3][M(HCOO)3] as Gua-M {i.e., Gua = 
guanidinium). 

Our general strategy for preparing [A][M(HCOO)3] sam¬ 
ples was as follows. Methanolic solutions of HCOOH ( 0.5 M, 
5 mL) and of the (usually) neutral A-site amine ( 0.5 M, 5 mL) 
were mixed at the bottom of a glass vial. Onto this solution 
methanol (2 mL) was carefully added, followed by a methano¬ 
lic solution of the transition metal nitrate ( 0.1 M, 8 mL).®The 
tube was then sealed and kept undisturbed. Following precipi¬ 
tation of the product (often as single crystals), the solution was 
filtered off, and the sample was washed with methanol, dried 
in air, and ground. The precise species used in the methanolic 
amine preparation were: methylamine (Acros Organics, 2 M 
solution in methanol), dimethylamine (Sigma Aldrich, 2 M so¬ 
lution in methanol), ethylamine (Sigma Aldrich, 2 M solution 
in methanol), azetidine (Sigma Aldrich, 98 %), and guanidine 
carbonate (Aldrich, 99 %). 
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Fig. 2 Compositions of the various ABX 3 formates used in our 
study. For the structural representations of the different A-site 
cations given at the top of the table, N atoms are shown in blue, C 
atoms in black, and H atoms in pink. The shaded regions of the table 
show the four compositional studies possible in which one or other 
of the A and B components is kept constant. 


2.2 X-ray powder diffraction 

Synchrotron X-ray powder diffraction data were collected us¬ 
ing the Ill beamlineES at the Diamond Light Source (A = 
0.82715 A) for each of the 19 different ABX 3 compositions 
given above. Finely-ground powder samples were loaded into 
0.5 mm diameter borosilicate capillaries and mounted on the 
diffractometer. The Cryostream Plus from Oxford Cryosys- 
tems was used to vary the temperature between 110 and 300 K, 
and diffraction patterns were collected using the Mythen2 po¬ 
sition sensitive detector (PSD). For each sample, data col¬ 
lection started at 300 K, before cooling to IlOK at a rate of 
5Kmin“^. To minimise the effects of beam damage, data 
were not collected continuously but rather at intervals of 10 K 
during this cooling process. Once the minimum tempera¬ 
ture was achieved, the goniometer was translated to allow a 
fresh part of the sample to be illuminated during heating. The 
sample temperature was then increased to 300 K at a rate of 
6 K min“\ with data collected at intervals of 10 K. For all data 
collected, diffraction patterns were obtained using two sepa¬ 
rate measurements (5 s each) at different angular orientations 
offset by 0.25° that were subsequently merged.l^Sl 

Structural models were determined from the powder 
diffraction patterns using Rietveld refinement as imp lemented 
in the TOPAS software (academic version 4.1).^^ For most 
of the systems studied, the corresponding crystal structures 
are already well known; the literature v alues wer e used as 
a starting model for Rietveld refinement.^^EHHEUE] 


best of our knowledge, the structures of MeNH 3 -M (M = Mg, 
Fe, Co, Zn, Cd) have not been determined previously; we 
found these to be isostructural to MeNH 3 -Mn and we used the 
known structure of this compound as a starting point for Ri¬ 
etveld refinement. In the case of Gua-Cd—which is also pre¬ 
viously uncharacterised—we made use of both single-crystal 
X-ray diffraction^ and our synchrotron X-ray powder diffrac¬ 
tion measurements to determine a relevant structural model 
(see SI for further details). 

Our Rietveld refinements made use of a number of 
chemically-informed restraints and constraints. Bonding re¬ 
straints were applied to the metal-formate bonds, while rigid 
bodies were used to model the organic molecular units (i.e. the 
A-site cations and the formate ligands). The bond lengths and 
internal angles defining the rigid bodies were refined for the 
initial temperature point. After refinement of the ambient- 
temperature structure for all [A][M(HCOO) 3 ] compounds, 
subsequent diffraction patterns were refined sequentially using 
the refinement model from the previous temperature as a start¬ 
ing model. In these sequential refinements, the free variables 
were: a polynomial background function, the lattice parame¬ 
ters, the atomic coordinates, a scale factor, and peak shape pa¬ 
rameters. The molecular geometries of the formate anions and 
the A-site cations were modelled as rigid bodies with transla¬ 
tional and rotational degrees of freedom. The Rietveld fits to 
the ambient and 110 K temperature points are given as SI. 


3 Results and Discussion 

3.1 Thermal expansivities and XBU analysis 

Our determination of the thermal expansion characteristics of 
ABX 3 formates is based on the interpretation of temperature- 
dependent lattice parameters extracted via Rietveld refinement 
of synchrotron X-ray powder diffraction data. While the var¬ 
ious phases we study are isostructural in the sense that their 
topologies are identical, they adopt a variety of different space 
group symmetries. The various crystal phases included in our 
study are shown in Fig.[^ What is immediately clear is that di¬ 
rect comparison of lattice expansivities from one system to an¬ 
other is not physically meaningful, since the various structures 
involve different relationships between the perovskite lattice 
and the unit cell geometry. 

Our solution to this problem is to interpret the lattice param¬ 
eter changes for a given system in terms of two fundamental 
XBUsiSIthe framework strut length r and the intra-framework 
framework angle 6. For systems of sufficiently high symme¬ 
try (e.g. the rhombohedral structures of Me 2 NH 2 -Mn/Co/Zn 
and Gua-Cd) there is a one-to-one mapping between the lat¬ 
tice parameters and XBU coordinates r, 0. In other words, the 
elastic behaviour of the framework is completely described by 
the propensity for hingeing (changes in 6) and network defor- 
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MeNH3-Mg EtNH3-Mn Me^NH-Mn Me^NH-Cu Aze-Mn Gua-Mn Gua-Cu Gua-Cd 

Pnma Pnl^a R3c CUc Pnma Pnna Pnl^a R3c 

Fig. 3 Relationship between unit cell and perovskite-like cubic structural unit for each of the crystallographically-distinct phases in our study. 
In each case, the transition-metal coordination polyhedra are represented as filled octahedra. C atoms are shown in black, N atoms in blue, O 
atoms in red, and H atoms in pink. 







mation (changes in r). For lower-symmetry structures, there 
will be more than one crystallographically-distinct value of r 
and/or 6 ; we use an average value in order to facilitate com¬ 
parison between different systems (see SI for further details). 
This means that the relationship between lattice parameters 
and r, 6 is approximate rather than exact, but we find that in 
practice the degree of approximation is < 10%, which is suf¬ 
ficient for our analysis. 

We demonstrate this approach using MeNHs-Mn as a rep¬ 
resentative example. This compound has the Pnma struc¬ 
ture shown in the leftmost panel of Fig. The relative 
changes in lattice parameters determined using Rietveld re¬ 
finement against our experimental X-ray powder diffraction 
are shown in Fig.|^ These reveal that, on heating, the frame¬ 
work expands rapidly along the ( 2 -axis and contracts almost as 
rapidly along the c-axis; this collective behaviour corresponds 
to hingeing of the framework as discussed elsewhere. ■3 The b- 
axis length is unaffected by framework hingeing, and its mod¬ 
est expansion with temperature reflects the intrinsic positive 
thermal expansion characteristic of M-formate-M linkages.® 
We quantify the magnitude of thermal response in terms of the 
coefficients of thermal expansion 


1 di 

ZdT’ 


( 1 ) 


to give aa = +88(3) MK \ = +19.5(4) MK ^ and = 

-49(2) MK-'P 

For this particular structure, the projection from lattice pa¬ 
rameter coordinates onto XBUs is given by the pair of equa¬ 
tions 


r = I + , 

0 = 2tan-'0. 


( 2 ) 

(3) 


The corresponding change in XBU magnitude with temper¬ 
ature is shown in the bottom panel of Fig. from which 



r(K) 


Fig. 4 Thermal expansion behaviour of MeNH3-Mn, represented in 
terms of the temperature dependence of lattice parameters (top 
panel) and XBUs (bottom panel). Data collected during cooling are 
shown as open circles; those collected during heating are shown as 
filled circles. The discrepancy between cooling and heating runs 
likely reflects a combination of beam damage and thermal lag. 
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Table 1 Lattice and XBU coefficients of thermal expansion for the 110-300 K temperature range, given in units of MK“^ For orthorhombic 
and hexagonal structures a\ = aa, cci = for monoclinic structures the at index represents the principal coefficients of thermal 

expansion.!^ The expansivities for Aze-Mn given here were determined for the orthorhombic cell characteristic of the ambient-temperature 
phase; see SI for the principal-axis expansivities of the low-temperature monoclinic phase. ^ For this phase, values of a were calculated using 
cooling data only. 


A+ 

M 2 + 

Space group 

ai 

a2 

a3 

ar 

oce 

MeNH 3 

Mg 

Pnma 

54.5(1.3) 

22.1(5) 

-20.0(5) 

20.1(5) 

46.0(1.1) 


Mn 

Pnma 

88(3) 

19.5(4) 

-49(2) 

21 . 8 ( 6 ) 

84(3) 


Ee 

Pnma 

74(2) 

14.7(4) 

-25.1(1.2) 

22.5(6) 

61(2) 


Co 

Pnma 

68.7(1.3) 

22.3(4) 

-28.8(6) 

21.6(4) 

61.1(1.2) 


Zn 

Pnma 

69(3) 

18.9(9) 

-34.6(1.0) 

18.9(1.2) 

65(3) 


Cd 

Pnma 

102(7) 

15.0(1.4) 

-61(7) 

21.5(7) 

100 ( 8 ) 

EtNH 3 

Mn^ 

Pn2\a 

44.2(4) 

33.2(4) 

4.14(13) 

28.4(3) 

23.88(15) 

Me 2 NH 2 

Mn 

R3c 

4.5(5) 

4.5(5) 

60.4(1.3) 

25.9(6) 

24.2(6) 


Co 

R3c 

10 . 6 ( 2 ) 

10 . 6 ( 2 ) 

46.5(1.0) 

24.2(5) 

15.5(3) 


Cu 

Cljc 

-14.3(1.0) 

45.7(1.7) 

57(3) 

31.7(1.6) 

24.7(6) 


Zn 

R3c 

8.3(9) 

8.3(9) 

36.4(4) 

18.9(1.9) 

12 . 2 ( 1 . 2 ) 

Aze 

Mn 

Pnma 

104(8) 

15.3(1.2) 

- 21 ( 6 ) 

31.6(1.1) 

-78(9) 

Gua 

Mn 

Pnna 

42.2(6) 

30.0(4) 

-10.6(3) 

19.2(3) 

-32.2(5) 


Ee 

Pnna 

33.1(4) 

31.6(3) 

-1-5(2) 

20 . 2 ( 2 ) 

-21.2(3) 


Co 

Pnna 

36.6(6) 

25.8(5) 

-6.7(2) 

17.5(3) 

-26.3(4) 


Ni 

Pnna 

27.1(8) 

22.4(6) 

3.59(16) 

17.1(5) 

-14.3(4) 


Cu 

Pnlia 

45.9(1.8) 

21.4(9) 

1.4(8) 

22 . 0 ( 1 . 1 ) 

-27.2(7) 


Zn 

Pnna 

30.7(1.9) 

21.9(9) 

-5.3(6) 

14.9(7) 

-21.9(1.4) 


Cd 

R3c 

-16.8(9) 

-16.8(9) 

106(3) 

16.3(3) 

-43.7(1.3) 


it is clear that the thermomechanical response is dominated 
by hingeing (changes in 6) rather than network deforma¬ 
tion. The XBU coefficients of thermal expansion are = 
+21.8(6)MK-' and ae = +84(3)MK"'. XBU projection 
equations for the various crystal structures of Fig. are given 
as SI. 

The extent of approximation in reducing the three lattice 
degrees of freedom to the two XBU degrees of freedom can 
be assessed by calculating values of a for the XBU-derived 
lattice parameters 



2 rtan( 0 / 2 ) 

(4) 

<^XBU 

V^l +tan2(0/2) 

^XBU 

= 2r 

(5) 


2r 

(6) 

<^XBU 

-yi +tan2(0/2) 

obtained via inversion 

of Eqs. (|^ and ([^: we find 

^«XBU — 

+87(3) MK-i, at^^^ 

, = +21.8(6) MK-i and 

^CxBV — 

—50(2) MK . These are all equal within error to the orig- 

inal lattice coefficients of thermal expansion. 



The ABX 3 lattice and XBU expansivities determined from 
our entire ensemble of X-ray diffraction data and correspond¬ 
ing Rietveld fits are summarised in Table the raw lattice pa¬ 
rameter data from which these values were derived are given 


as SI. While the magnitudes of the lattice expansivities vary 
over two orders of magnitude, we find essentially the same ba¬ 
sic thermomechanical response for all ABX 3 structures. For 
example, nearly all of the orthorhombic structures exhibit 
large positive thermal expansion (PTE) along the ( 2 -axis (27- 
104MK“^), moderate PTE along the b-nxis (14.7-33 MK“^), 
and NTE along the c-axis (—61 to —1.5 MK“^). 


3.2 Effect of variation in 

In order to investigate the effect of variation on the me¬ 
chanics of [A][M(HC 00 ) 3 ] frameworks, we compare our re¬ 
sults for the MeNH 3 -M and Gua-M families—these are the 
two systems for which we have the greatest diversity of 
substitution. Eig. [^presents graphically the relevant data from 
Tabl^ organised according to the (Shannon) radii of the 
ions.^Two results are immediately obvious from considera¬ 
tion of the dependence of XBU expansivity on ionic radius. 
The first is that the value of is essentially composition- 
independent; the second is that the magnitude of ae increases 
linearly with increasing radius. 

The magnitude of reflects the effect of increased ther¬ 
mal motion on the M.. .M separation across connected M- 
formate-M links. Our Rietveld refinements do not allow us 
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Lattice XBUs 




Fig. 5 Lattice expansivities as a function of metal cation size for (a) MeNH3-M and (b) Gua-M families (only data for the orthorhombic 
members are shown). The dotted ellipse in (b) indicates the aa value for Gua-Cu, which is not included in the linear fit shown here as a 
guide-to-the-eye. XBU expansivities calculated for complete (c) MeNH3-M and (d) Gua-M families. Ionic radii are taken from Ref. [ 57 lfor 
octahedrally coordinated (high-spin) cations. 
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to comment with certainty on the temperature-dependence of 
the position, orientation, and/or thermal displacement of the 
formate anions in the various compounds in the way that 
has allowed detailed analysis of thermal expansion effects 
in other dense MOFs.^ Nevertheless a likely origin of the 
composition-independence of is the balance between two 
competing effects: on the one hand, larger ions would 
allow greater vibrational motion of bridging formates (i.e., 
lower vibrational frequencies); on the other hand, the in¬ 
creased M.. .M separation for larger ions means that formate 
displacements have less effect on this separation in relative 
terms. So while it is likely that the thermal motion of formate 
ions increases with increasing radius, the relative effect 
on r as an XBU remains roughly constant. 


By contrast, the value of represents the degree of hinge- 
ing flexibility in a given framework. The increase in flexibil¬ 
ity we observe with increasing cation size could arise from 
several effects. First, the strength of the metal-formate coor¬ 
dination bond will decrease as larger metal cations are used, 
which in turn facilitates the movement of this bond. This 
causes an enhancement of the framework flexibility towards 
framework hingeing, especially as the neighbouring formate 
linkers are now further apart within the coordination sphere of 
the metal cation. In previous compositional dependent stud¬ 
ies of Prussian Blue analogues, the strength of the coordina¬ 
tion bonds was important for tuning the transverse vibrations 
of the M-ligand-M strut.® Second, the larger unit cells ob¬ 
tained with larger metal cations may allow greater structural 
freedom, especially as the hydrogen bonding from the A-site 
cations to the anionic formates of the framework are at greater 
distances. These two factors mean that small metal cations, 
such as Ni^+, will have (/) short and strong M-0 bonds, re¬ 
stricting M-formate movement, and (//) a small unit cell, thus 
forming shorter intermolecular contacts between formate link¬ 
ers, and between the framework itself and the A-site cation; all 
of these factors will result in a restriction of framework hinge¬ 
ing. 


While the values of are essentially identical for both 
MeNHs-M and Gua-M families 20MK“^), the magni¬ 
tudes of aQ and their sensitivity to cation radius are substan¬ 
tially smaller for the latter family than for the former. We dis¬ 
cuss this dependence on alkylammonium cation in more detail 
below, but note here that this difference has the consequence 
for Gua-M systems with small that | | ^ | | and so the 

NTE effect that would ordinarily arise from framework hinge¬ 
ing can be masked by the network deformation (expansion in 
r). This is why Gua-Ni and Gua-Cu do not exhibit NTE along 
any crystal axis. 


Table 2 Computed principal axis effective lengths (given in units of 
A^) for each of the A-site cation centr oids i n the host manganese 
formate, calculated within CRYSTALS.®^ The asphericity h (also 
given in A^) is defined by Eq. 0. 


A+ 

Lmin 

f^med 

f^max 

b 

MeNH 3 

0.54 

0.57 

1.88 

-0.64 

EtNH3 

0.65 

0.88 

2.95 

-0.92 

Me 2 NH 2 

0.71 

1.21 

3.25 

-0.77 

Aze 

0.82 

1.92 

1.98 

0.52 

Gua 

0.48 

2.39 

2.52 

0.89 


3.3 Effect of variation in A+ 

Our results for the MeNH 3 -M and Gua-M families show 
clearly that variation in A-site cation can affect the extent of 
framework hingeing observed. To investigate this in more 
detail, we consider the family of Mn-containing ABX 3 for¬ 
mates with A = CH3NH3+, CH3CH2NH3+, (CH3)2NH2+, 
(CH 2 ) 3 NH 2 +, and C(NH 2 ) 3 +. The sizes of the different A- 
site cations were estimated by fitting the atomic coordinates 
refined from the powder diffraction data to a shape tensor, L, 
using the program CRYSTALS. ® This shape tensor represents 
an anisotropic ellipsoid that contains all atoms within the A- 
site cation.® The smallest, median, and largest components of 
this ellipsoid are shown in Table We initially use the maxi¬ 
mum effective length of the A-site cation (Lmax) as our metric 
by which to compare the mechanical properties of our differ¬ 
ent systems. In a another study, a different method was used to 
calculate the sizes of the A-site cations mentioned here.® The 
two approaches are consistent with one another except that the 
order of the Me 2 NH 2 ^ and Gua+ cations is reversed (see SI 
for further discussion). 

Figure shows the XBU expansivities for the five Mn- 
containing formate frameworks, arranged according to the 
value of Lmax for their A-site cation. As observed in the previ¬ 
ous section, there is essentially no meaningful compositional 
dependence of a^. By contrast, the magnitude of is much 
reduced for those systems with larger A-site cations. Our use 
of the absolute value of reflects the fact that its sign is not 
conserved amongst the various compounds we study, a point 
we will expand upon below. So, at face value our results sug¬ 
gest that bulkier A-site cations inhibit framework flexibility. 
This may be due simply to steric interactions between the A- 
site cation and the host framework, with bulkier cations pre¬ 
venting large changes in framework angles. From a chemical 
viewpoint it is likely that the different H-bonding strengths of 
the various A-site cations will also impact the framework flex¬ 
ibility;® however, our data do not allow us to probe changes 
in hydrogen-bond characteristics with any certainty. Even in 
these simplistic terms, and taken together with the results of 
the previous section, our analysis already suggests that the 
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Fig. 6 XBU expansivities determined for MeNH3-Mn, Aze-Mn, 
Gua-Mn, EtNH3-Mn, and Me2NH2-Mn, given as a function of 
maximum effective A-site cation length. The cation shape tensor 
generated using CRYSTALS is represented for each compound below 
the graph.l^^ES The expansivities for Aze-Mn have an anomalously 
large uncertainty as a result of a low-temperature phase transition 
(see SI for further discussion). 


most flexible formate frameworks will be those with large 
but small A+—precisely the combination that (rightly) 
identifles MeNH 3 -Cd as the most flexible of the systems we 
study here [Table 

One of the key motivations for our study was to under¬ 
stand the relationship between A-cation shape and framework 
mechanics. There are a large number of metrics to quantify 
shape, but we constrain ourselves here (given that we have but 
flve data points to compare) to the straightforward notion of 
asphericity^T\YQ asphericity parameter 
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Fig. 7 Relationship between asphericity b and inverse XBU 
expansivity for the five Mn-containing ABX3 formates in our study. 
The direction of framework flexing switches for prolate/oblate 
A- site cations; the extent of fiexibility is largest for small, spherical 
cations (Z? ^ 0). 


ae > 0 and ao < 0, respectively!^ This switching in sign of 
anisotropic mechanical response as a consequence of cation 
asphericity is strongly reminiscent of the critical ratios that de¬ 
marcate linear and area negative thermal expansion in a range 
of molecular framework materials.®! As in those systems, the 
magnitude of hingeing response is actually larger for systems 
closer to the critical geometry (here, /? = 0). Our data show an 
approximately-linear relationship between b and the inverse 
expansivity ^ [Fig. . 

So, from a design perspective, our analysis suggests that it 
is not only the A-cation size that affects framework flexibility, 
but also its anisotropy. Hence we can explain why EtNH 3 -Mn 
is less flexible than Aze-Mn, despite the fact that EtNH^ is the 
smaller cation by ^25% (see SI for details). This relationship 
also explains why the qualitative ordering of Fig.[^in terms of 
^max is physically meaningful. 

4 Concluding Remarks 


^ — ^med 2 (^min ^max) 


( 7 ) 


takes into account both the size of an object and the anisotropy 
of the shape tensor L. The parameter b assumes negative val¬ 
ues for objects with prolate asphericity, positive values for 
oblate objects, and is zero for isotropic shapes. The var¬ 
ious A-site cations we consider span a range of aspheric- 
ities — 1 ^ b < 1A^ [Table i. Re-ordering the hingeing 
XBU expansivities according to these values, we And a clear 
distinction between prolate and oblate A-cations, for which 


The relationships we establish between cation size and shape 
provide a set of straightforward structure-property rules 
which can be implemented as part of the strategic design of 
functional MOFs. The basis of these relationships is the cor¬ 
relation between cation size and M-formate bonding strength. 


* Note that in order for the sign of aQ to have any physical meaning, we are 
careful to define the value of 0 in a consistent way for all frameworks we 
study. In particular, it is given as the obtuse angle in the rhombic face of the 
perovskite-like cube. Hence > 0 implies that the cube geometry becomes 
more distorted with increasing temperature, and < 0 reflects a tendency to 
adopt an increasingly cubic geometry on heating. 
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which in turn affects the mechanical properties. In a com¬ 
putational study of the [Me 2 NH 2 ][M(HCOO) 3 ] family, it was 
shown that the calculated Bader charges on the metal cations 
exhibited the following order: Mn > Fe > Zn > Co > Ni.l^ 
This order was also suggested to represent the M-formate 
bonding strength, where Mn-formate exhibits the weakest 
bond (highest ionic component) and Ni-formate the strongest 
(greatest covalent component), and furthermore, the ordering 
is very similar to that of the metal cation size. There is a small 
discrepancy with the Zn and Co placement between the Bader 
charge and metal cation size sequences; however, as there is 
very little difference in the sizes (0.74 0.745 A), the varia¬ 

tion is not significant. This small discrepancy can be also seen 
in the coefficients of thermal expansion, which vary subtly in 
the order for these two metal cations [Fig.[^. 

In a nanoindentation study of the [Me 2 NH 2 ][M(HCOO) 3 ] 
family, it was suggested that the metal ligand field stabilisation 
energy (LFSE) directs the resulting mechanical properties. E2 
The LFSE order of Mn = Zn < Co < Ni was reproduced in 
the magnitudes of the Young’s Moduli of these compounds.!^ 
This finding places compounds which contain either Mn^+, 
Zn^+, or Cd^+ on the same level of structural fiexibility, which 
is not consistent with the thermal expansivity characteristics 
we measure here. Likewise, a Brillouin scattering study of 
related [NH 4 ][Mn(HCOO) 3 ] and [NH 4 ][Zn(HCOO) 3 ] struc¬ 
tures indicated increased stiffness for M = Zn^+ than M = 
Mn^+PThe thermal expansion data from this study and oth¬ 
ers® also show increased stiffness for Zn^+-containing com¬ 
pounds compared to Mn^+. This discrepancy could be due 
to the different experimental methods used: in the case of 
nanoindentation experiments, a uniaxial force is exerted upon 
the crystal, while upon temperature variation, an isotropic ex¬ 
ternal stimulus is applied. One obvious and straightforward 
experiment would be to carry out a comparative nanoinden¬ 
tation study of Zn- and Cd-containing formate frameworks, 
since these species correspond to very different ionic radii but 
identical LFSEs. 

A clear result of our study has been to show that the pres¬ 
ence of A-site cations within the pores of the metal formate 
framework plays an important role in controlling the struc¬ 
tural fiexibility of the material. In particular, longer A-site 
cations cause a significant decrease in the framework hinge- 
ing observed. Thus these structures are not expected to give 
rise to anomalous mechanics, such as NTE. Instead their me¬ 
chanics will be dominated by the behaviour of the M-linker- 
M units. In addition, the direction of framework hingeing can 
be switched by using differently-shaped A-site cations within 
the framework pores. In the case of oblate A-site cations, 
the framework hingeing is directed towards a convergence of 
framework angles (such as 90° for a 2D wine-rack) upon heat¬ 
ing; whereas there is a divergence of framework angles in 
the case where prolate A-site cations are used. The frame¬ 


work hingeing direction could have implications for ferroelec¬ 
tric transitions arising from reorientation of extra-framework 
molecular ions. Stronger host-guest interactions might be ex¬ 
pected for pore shapes which mimic the shape of the molecu¬ 
lar ion. Thus ferroelectric ordering may be stabilised in cases 
where the pores vary (Le. framework angle evolve) towards a 
shape that mimics that of the molecular ion. 

In summary, our work has shown that metal cation size cor¬ 
relates well with expansivity magnitudes. Eor the octahedral 
cation site, cations with smaller radii give rise to stiffer frame¬ 
works, while the larger metal cations are associated with the 
more extreme framework fiexibility. At the A-site, cation size 
is also found to affect framework fiexibility, though in the op¬ 
posite sense: framework hingeing is constrained when long 
A-site cations are present within the framework pore. More¬ 
over, the type of asphericity of the A-site cation determines 
the direction of framework hingeing. This is most likely to be 
observed only in dense MOEs, where the vibrational motion 
of the anisotropic A-site cation is of greater importance. The 
most readily applicable rule for rational mechanical design of 
MOEs involves that of the metal cations, where simple evalua¬ 
tion of its size within different oxidation states or coordination 
environments can provide a scale of fiexibility within a series 
of isostructural frameworks. 
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* Single-crystal X-ray diffraction data for Gua-Cd were collected using a 
Nonius KappaCCD diffractometer or an Oxford Diffraction (Agilent) Su- 
perNova diffractometer fitted with an Oxford Cryosystems Cryostream 600 
Series/700 Plus open flow nitrogen cooling device. denzo/scalepaciJ^ 
or CrysAlisPro were used for data collection and reduction as appropri¬ 
ate. In general, the structures solved ab initio using SIR92I^ although co¬ 
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tures. All stru ctures were refined with full-matrix least-squares on us¬ 
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ence Fourier map and treated in the usual manner.!^ Full structural data 
are included in the Supplementary Information (SI) and have been submit¬ 
ted to the CCDC as numbers 1420162-1420165. These data can also be ob¬ 
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